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Asymptotic similarity states at large Reynolds numbers and small Rossby num- 
bers in rotating homogeneous turbulence are investigated using the database ob- 
tained from large-eddy simulations of the incompressible Navier-Stokes equations. 
Previous work has shown that the turbulence kinetic energy and integral length 
scales are accurately described by simple scaling laws based on the low wavenumber 
part of the three-dimensional energy spectrum (Squires et al. 1994). The primary 
interest of the present study is to search for spectrum similarity in the asymptotic 
state. Four independent energy spectra axe defined. It is shown that rescaling of 
these energy spectra in the asymptotic regime will collapse three out of the four 
spectra. The spectrum which does not collapse is a function only of the vertical 
wavenumber and corresponds to two-component motions in the plane normal to 
the rotation axis. Detailed investigation of the cause of this anomalous behavior 
reveals the existence of a strong reverse cascade of energy from small-to-large scales 
of the two-dimensional, two-component motions. This feature of the rotating flow 
is presumably linked to the lack of a complete similarity state, though further study 
of this issue is required. 


1. Introduction 

Solid-body rotation of initially isotropic turbulence represents the most basic tur- 
bulent flow whose structure is altered by system rotation but without the complicat- 
ing effects introduced by mean strains or flow inhomogeneities. Squires, Chasnov, 
Mansour & Cambon (1994) (referred to as SCMC throughout this manuscript) ex- 
amined the long-time, asymptotic evolution of rotating homogeneous turbulence at 
large Reynolds numbers and small Rossby numbers. They found that the kinetic 
energy decay exponent was reduced by approximately a factor of two relative to its 
value in non-rotating turbulence and was independent of the rotation rate. SCMC 
also found that the integral length scales undergo power-law growth and that the 
evolution of both the kinetic energy and length scales in the asymptotic regime is 
accurately described using simple seeding laws based on the invariants of the low 
wavenumber part of the initial energy spectrum. 
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The scaling laws of the kinetic energy and integral length scales deduced by 
SCMC suggest the existence of an asymptotic similarity state in rotating homoge- 
neous turbulence. Since a similarity state permits prediction of the flow evolution 
using a simple rescaling of the statistics, existence of asymptotic similarity states 
enables prediction of the ultimate statistical evolution of the flow without detailed 
knowledge of non-linear transfer processes. Similarity states thus imply a certain 
universal character to a flow which may ultimately simplify modeling of complex 
turbulent flows at high Reynolds numbers. 

The objective of the present work, therefore, is to search an asymptotic similarity 
state for the energy spectra. In §2 an overview of the LES database used in this 
study is summarized. Also presented in §2 is the time development of the kinetic 
energy power-law exponent, integral length scales, and asymptotic scaling laws ob- 
tained from both dimensional analysis and the simulation database. The reader is 
referred to Squires et al. (1993, 1994) for a more complete discussion. In §3 the 
asymptotic spectrum similarity is analyzed and tested using the LES database. 

2. Database overview and scaling laws 

2.1 Database overview 

The database used in the present study was obtained using large-eddy simulation 
of the incompressible Navier-Stokes equations in a rotating frame of reference. The 
use of a large-eddy simulation with zero molecular viscosity ensures that the flows 
considered are nominally at infinite Reynolds numbers. For purposes of discussion 
the vertical (r 3 ) axis is taken along the axis of rotation. The governing equa- 
tions were solved using the pseudo-spectral method developed by Rogallo (1981). 
Subgrid-scale motions were parameterized using a spectral-eddy viscosity modified 
for system rotation (Kraichnan 1976, Chollet & Lesieur 1981, Aupoix 1984). Cal- 
culations were performed from isotropic Gaussian initial conditions using an initial 
energy spectrum of the form 


E{ Jfc,0) 


where C a is given by 



( 1 ) 

( 2 ) 


and k p is the wavenumber at which the initial energy spectrum is maximum. SCMC 
performed simulations with s = 2 and s = 4, corresponding to the initial energy 
spectrum having a low wavenumber form proportional to either k 2 or k 4 . 

Simulations were performed using resolutions of 128 x 128 x 512 on a computa- 
tional domain which is four times longer along the rotation axis than in the other 
directions. This yields an isotropic grip in physical space, but permits the de- 
velopment of large-scale anisotropies without compromising the usefulness of the 
computational results. The maximum wavenumber in each of the three coordinate 
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Figure 1. Time development of the power-law exponent of (u 2 ) in rotating 

turbulence, A; 2 spectrum: o © Q = 0; fi = 0.5; f2 = 1.0. k spectrum. 

ft = 0; ft = 0.5; n = 1.0. 


directions was 95, the initial root-mean-square velocity fluctuation u 0 was 1/2, and 
the wavenumber at which the initial spectrum was initially maximum, kp, was 75. 
For each spectrum type, i.e., low wavenumber part proportional to A; 2 or A: 4 , sim- 
ulations results are available for three rotation rates, fi = 0, 0.5, and 1.0. The 
flows were evolved over a sufficiently long period so that the effects of rotation on 
the flow dynamics were clearly evident, and the kinetic energy decay and integral 
length- scale growth approximated power-law forms. This corresponds to an evolu- 
tion of the flow to low Rossby numbers, where the Rossby number is a measure of 
the relative strength of inertial forces to Coriolis forces. 

2.2 Statistical evolution and scaling laws 

The effect of rotation on the evolution of the turbulence kinetic energy ((u 2 )) can 
be illustrated through examination of the power-law exponent, n (in (u 2 ) oc t n , 
see Fig. 1). In Fig. 1, the time axis is made dimensionless using the eddy turnover 
time in the initial field 

r( 0) = X u (0)/{u 2 ) 1 / 2 (3) 

where L u {t) is the velocity integral scale at time t defined as 


7T / 0 °° k~ l E(k,t)dk 
~2 J 0 °°E(k,t)dk 


( 4 ) 


Power-law exponents have been plotted for each initial spectrum type and rotation 
rate used in the calculations. It is evident from Fig. 1 that following an initial 
transient the power-law exponent becomes independent of time. Comparison of the 
power-law exponents for both initial spectrum types demonstrates the reduction 
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in the decay rate of kinetic energy in rotating turbulence. It is also clear from 
Fig. 1 that in the asymptotic regime the power-law exponent of the kinetic energy 
is independent of the rotation rate provided the rotation is non-zero. The value of 
the exponent depends only on the form of the initial energy spectrum. 

It is possible to derive the appropriate asymptotic scalings of (u 2 ) for non- rotating 
turbulence using dimensional analysis (e.g., see Chasnov 1994 ). For the rotating 
flow, assuming that non-linear triadic interactions are correlated over a time scale 
directly proportional to ft, SCMC showed that the asymptotic scalings for the 
kinetic energy are 

(u 2 ) ex I?Jr"s(ftf) T (fc 2 spectrum) (5) 

(u 2 ) a B% (k 4 spectrum) . (6) 

For the k 2 spectrum B 0 is invariant while for the k 4 initial conditions B 2 is time 
dependent. However, the time dependence of B 2 is weak relative to the overall 
turbulence decay. It is assumed invariant for the purposes of developing asymptotic 
scalings. The scaling laws ( 5 ) and (6) predict that in rotating turbulence the kinetic 
energy decay exponent is reduced by a factor of two, relative to its value in the non- 
rotating flow. As shown by Fig. 1, simulation results agree with these predictions. 

SCMC also examined the development of the integral length scales. These length 
scales are obtained from integration of the two-point correlation function 

Lacrf = j (u a (x)u a (x + rep))drp , (7) 

where e$ is the unit vector in the (5 direction, and a, j3 = 1, 2 , or 3 . The length scale 

L aa j measures the correlation between the a velocity components with separation 
in the 3 direction. SCMC examined the evolution of five length scales. In the 
horizontal plane these integral scales are 

Lhi = — (X-11,1 + £22,2), L h2 = ^(£11,2 + £22,1), L h3 = ^(£33,1 + £33,2). ( 8 ) 

The integral scale Lhi measures the horizontal correlation of longitudinal horizontal 
velocities, Zj, 2 the horizontal correlation of lateral horizontal velocities, and Lh 3 the 
horizontal correlation of vertical velocities. For vertical correlations (along the axis 
of rotation) there Eire two independent integral scales 

L v i = -(£11,3 + £22,3)1 L V 2 = £33,3 , ( 9 ) 

where L v % measures the vertical correlation of horizontal velocities and L v 2 the 
vertical correlation of vertical velocities. 

Shown in Fig. 2 is the time development of the integral scales from SCMC for the 
highest rotation rate used with each initial condition type. Also shown for reference 
is the slope corresponding to the asymptotic growth rate in non-rotating isotropic 
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FIGURE 2. Time development of integral length scales in rotating turbulence, 
fl = 1. Asymptotic growth rate in non-rotating turbulence also shown in lower 

portion of figure. Lhi) Lk2] ^vi] ( a ) 

spectrum; ( b ) k A spectrum. 


turbulence. It is clear from Fig. 2 that L vX grows much more rapidly than all other 
integral scales. Also, the growth of the other four integral scales appear to be slightly 
suppressed relative to the non-rotating case. In addition, the simulations provide 
evidence that all the integral length scales exhibit power-law growth asymptotically 
and, with the important exception of L vX , the power-law exponents are all the same. 
Thus, in the horizontal plane there appears to be a single independent length scale, 
which we denote generically as Lh, whereas there exists two independent length 
scales in the vertical direction, namely L v \ and L v 2 . SCMC further demonstrated 
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that the evolution of the integral scales from two simulations of a k ^ spectrum at 
different rotation rates fi could be collapsed according to the asymptotic scaling 
laws 

L v2 oc (fit)'* , (10) 

while for the A; 4 spectrum a collapse was achieved using the seeding laws 

z*oc* a M(nt)“*, l v1 l v2 * B*t* (at)'* . (ii) 

3. Asymptotic spectrum similarity 

Power-law behavior of both the kinetic energy (Eqs. 5 and 6) and integral length 
scales (Eqs. 10 and 11) suggest the existence of energy spectra which are asymptot- 
ically similar, i.e., the time- and rotation-dependent spectra collapse by appropri- 
ately rescaling the wavenumber and spectral density axes. For brevity, we consider 
here only our results for the k ^ initial conditions. More complete results will be 
presented in a future publication. 

As the first obvious candidate for an asymptotically similar spectrum, we define 
an axisymmetric energy spectrum as a function of k± and Ar 3 as 

$(*±>* 3 ,*) = ( 12 ) 

where 

k ± = (&? + k l)* (13) 

is the horizontal wavenumber perpendicular to the rotation axis, t/ * (k) is the Fourier 
transform of a periodic velocity field u,(x), * denotes the complex conjugate, and 
the angular brackets used in (12) denote an average over a horizontal ring in three- 
dimensional wavespace with center located at a fixed value of jfc 3 . The relationship 
between the spectrum in (12) and the kinetic energy per unit mass of the fluid is 
given by 

\( u2 ) = f dk * J q dk ± <Sf(k ± ,h,t). (14) 

An attempt to construct a similarity form directly for the spectrum <I>(fcj_, k 3 ,t) 
is impeded by the existence of two distinct vertical integral scales in the rotating 
flow. The choice of an appropriate length scale for rescaling k 3 is thus ambiguous, 
making it unlikely that the full spectrum ^(k±,k 3 ,t) will be asymptotically similar. 
We may, however, integrate out the vertical wavenumber from the spectrum defined 
in (12) to construct an energy spectrum which is a function only of the horizontal 
wavenumber: 

/ OO 

dk 3 $(k ± ,k 3 ,t). (15) 

-OO 

The simulation results have shown that there is only a single horizontal length scale 
(10) which may be used to scale the horizontal wavenumber This allows us to 
propose the following general form for the similarity state of the spectr um 'kjj 

*j.(k±,t)=(u 2 )L h * ± (L k k x ). 


(16) 
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FIGURE 3. Time development of ft = 0.5; ft — 1.0. (a) 

unsealed; (b) scaled. 


Making use of the known asymptotic behavior of the kinetic energy (5) and integral 
scales (10), the similarity state (16) is given explicitly as 

*x(kut)=B$t-t(Qt)i$±.(kx), k ± = . (17) 

Shown in Fig. 3 are the unsealed and scaled spectra 
after the flows have evolved sufficiently to be considered asymptotic in the sense that 
rotation effects are important (low Rossby numbers) and the statistics follow power- 
law forms. Spectra from calculations using both non-zero rotation rates have been 
included in the figures. While the decay with increasing time and the dependence 
on rotation rate is evident in Fig. 3a, the scaled spectra and wavenumbers (17) 
in Fig. 3b show a remarkable collapse, indicating that ^x is independent of both 
time and rotation rate. The results shown in Fig. 3b thus confirm the similarity 
state given by (17). This is a significant result since we have succeeded in reducing 
the number of independent variables in the problem from four, (&x> t, ft* and B o), 
where Bo represents the effect of initial conditions, to simply one, (k±). 

We must still confront the difficulty of finding similarity spectra which are func- 
tions of the vertical wavenumber k^. We know from the LES results that two 
vertical length scales exist, one associated with the horizontal velocities and the 
other with the vertical velocities. This suggests splitting the energy spectrum into 
two spectra, each separately related to the horizontal or vertical velocities. A sim- 
ple splitting, however, does not take into account the constraint imposed by the 
continuity equation on the flow field. Rather, we proceed to define spectra (Herring 
1974) by first projecting the Fourier components of the velocity field onto two unit 
vectors perpendicular to k (the velocity component parallel to k is zero by virtue 
of the continuity equation): 


U> (k) = ^ 1 (k)e! 1) (k) + ^(k)e; 2) (k), 


( 18 ) 
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where 


e (1) (k) 


k x n e W (k) _ k x e(1) ( k ) 

|k X nr W - |k X e0)(k)| • 


( 19 ) 


We note that the field <^i(k) consists only of horizontal velocity fluctuations, whereas 
the field <fo(k) contains both horizontal and vertical fluctuations. We then use the 
scalar fields 0i (k) and 02 (k) to define energy spectra in analogy to (12): 


* (1) (*j.,* 3 ,t) = (20) 

* (2) (*l, h,t) = k, t)<t> 2 { k, ty), (21) 

where we have now separated the energy spectrum into two spectra: 


*(*L,*3,<) = ¥ (1) (*x,*3,<) + ¥ (2) (*l,*3,0- 


( 22 ) 


The significance of this separation can be seen by expressing the integral length 
scales L vX and L v2 of (9) in terms of and i.e., 


7T f°° 

dkx¥ "^' 0) ' < 23 ) 

and 

= (24) 

The LES results show that the horizontal and vertical mean-square velocity fluctu- 
ations decay asymptotically with the same power-law form, so that (23) and (24) 
serve the hopeful purpose of matching the two independent vertical integral scales 
to the two different spectra. 

In analogy to the similarity form given in (16), we now postulate the following 
general similarity states of the form 


* {1) (k u k 3 ,t) = (25 a) 

* (2) (*.l,*3,0 = (<Pl)l h 2 l V 2* (2) (l%k ± ,l v 2 k 3 ), (25 b) 

where ($), /f, q, and /*, q are the mean-square fluctuations and horizontal 
and vertical integral scales associated with the ^ (1) and spectra, respectively. 
The relations between the mean-square fluctuations, the integral scales, and the 
spectra are given for reference by 

<tf,2> = I" dk 3 /°°d*x¥ (, ’ 2) (*L,*3,<), (26) 

J — co Jo 

and 


/ i 0 ’ 2 "< 4 , 2 )L c dk3 l 


* il ' 2) (kj_,k 3 ,t) 


( 27 ) 
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We have computed the amplitudes and integral scales given by (26) and (27) from 
the LES database. The mean-square values ( <f >\ ) and {4>\) decay asymptotically 
at the same rate as the total kinetic energy per unit mass as given in (5), and the 
horizontal integral scales Zj and l 2 and vertical integral scales Z 2 and 1% grow asymp- 
totically as Lh (for both horizontal integral scales) and L v \ and L v 2 , respectively, 
as given by (10). 

Although the similarity states proposed in (25) appear to be reasonable, a closer 
examination reveals a further problem. We recall that the asymptotic form of the 
energy and integral scales were constructed on the basis of the invariant Bq , related 
to the form of the energy spectrum near zero wavenumber. Assuming a flow field 
starting from isotropic initial conditions (as in our simulations), the forms of the 
spectra ^ (1) (^’x, h,t) and <P (2) (fcx, h,t) near k = 0 are 


* (1 \kx,k 3 ,t) = 2nk ± (B 0 + O(k)), 

* w (k±,k 3 ,t) = 2*k ± {B 0 +O(k)), (28) 

1/2 

where O(Jfc) represents terms which vanish as k = + fc 2 ) — * 0 and Bo is 

independent of time. Near k = 0 the similarity forms of the spectra must necessarily 

bC $ w (l h l k ± ,l v l kz)~2nl h 1 k ± , $ (2) (Z 2 Hx,/ 2 0 M~27rZ 2 fc Ar x ; (29) 

and substituting (28) and (29) into (25) yields 

Bo oc (</>?)(Zl*) 2 Zr , Bo oc (^)(Z 2 fc ) 2 / 2 ° • (30) 

Eq. (30) is an important result in that it provides two additional constraints on 
the allowed form of an asymptotic similarity state based on the invariant B 0 . The 
second constraint of (30) can be seen to be satisfied by the LES results, where 

z 2 h ocfl 0 M(m)-*, q <xB 0 M («*)"* . (3i) 

However, the first constraint of (30) is widely violated by the LES results due to 
the rapid growth in the vertical integral scale: 

(^ocBoh-qn*) 1 , /f<x* 0 M(nt)-*, ijocflftqn*) 1 •> ( 32 ) 

so that 4 

(4> 2 M) 2 n cxBo(ta)*. 

The similarity state in (25a) constructed for t) must thus be in error. 

However, the similarity state in (25b) constructed for *) ma y possibly 

be correct in so far as the constraint given by the second of (30) is satisfied. 

To test the correctness of the similarity state for we ^ s °l a ^ e the 

troublesome jfc 3 dependence of the spectrum by integrating over k± to obtain the 
(usual) one- dimensional spectrum 
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Figure 4. Time development of 4>3 2) (& 3 ). ft = 0.5; ft = 1.0. (a) 

unsealed; (6) scaled. 



Figure 5. Time development of ft = 0.5; ft = 1.0. (a) 

unsealed; ( b ) seeded. 


*r(* 3 ,<)= / dk ± ¥ 2 \k u h,t)-, (33) 

Jo 

and the corresponding similarity state is postulated to be 

= (41)1^1(1^), (34) 

which using (31), may be written explicitly as 

*s 2) (*»»0 = B 0 V«(ftt)«$< 2) (£ 3 ), h = 3 . (35) 

Shown in Fig. 4 are the unsealed spectra ’^(A*,*) and the scaled spectra 
$< 2) (£ 3 ,f) obtained from (35). As before, results for both non-zero rotation rates 
are shown in the figure. In Fig. 4 a the spectra at different times and rotation rates 
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FIGURE 6. Time development of £3) (Fig. 6a) and ,£3) (Fig 6b) 

for A?3 — 0 and fl = 1.0. Increasing time is from to e » ; asymptotic 

behavior as k± — » 0 (Eqn. 29). 


are clearly distinguishable whereas in Fig. 46 the spectra for the times 

plotted and both Q collapse onto a single curve. Together with the previous results 
of Fig. 3, we have thus confirmed the complete similarity state of ^ 2 ^(fc_L, ^3, t) 
given by (25b). 

For comparison, the troublesome one-dimensional spectrum and the 
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scaled spectrum corresponding to the assumed similarity form given by (25a) are 
shown in Fig. 5. As anticipated from the inconsistency of the assumed similarity 
form with the constraint given by the first of (30), the rescaling shown in Fig. 56 
does not collapse 4/ 3^ (£3) onto a single curve. 

Thus we have demonstrated spectral similarity for three out of four spectra, 
namely and are asymptotically similar, whereas 

is decidedly not. 

Further evidence of the anomalous asymptotic behavior of ^3^ is provided in 
Fig. 6. Shown in the figure, for Q = 1.0, are tyW (Fig. 6a) and 4d 2 ) (Fig. 66) 
as a function of k± for = 0. From (28), the spectra in these figures should 
follow the asymptotic form \t (1)2 >(fci, ir 3 - 0 ,*) ~ 2nk±B Q as -> 0. Considering 
Fig. 66, 0, /) decays as one would expect from a decaying three-dimensional 

turbulence. The energy at larger values of k± decay more rapidly than the energy 
at smaller values, and consequently the peak of the spectrum systematically moves 
to lower wavenumbers. The invariance of the form of the spectrum as k± — 1 ► 0 is 
respected. The spectrum 0,t) shown in Fig. 6a, however, behaves much 

differently. A strong increase in the energy at small values of k± is clearly evident. 
The invariance of the spectrum as k± — ► 0 is a rigorous result, yet appears to 
be irrelevant to the time- development of this spectrum. It is important to note 
that the spectrum t) corresponds to two-dimensional, two-component 

motions, and what we are observing in Fig. 6a is apparently a reverse cascade of 
energy from small to large scales, well-known within the context of two-dimensional 
turbulence. The existence of a reverse cascade may be the reason why a similarity 
state based on the invariance of the low wavenumber coefficient does not develop 
for the spectrum. The establishment of a two-dimensional reverse cascade 
is an important feature of three-dimensional turbulence in a rotating frame, and 
obviously requires a deeper understanding than we have achieved thus far. Other 
recent investigations (Bartello ei al 1994; Hossain 1994) have also confronted this 
very question, although from a different perspective than the present work. 

Finally, the strong reverse cascade of energy of the two-dimensional, two-compo- 
nent motions has interesting implications as to the asymptotic structural evolution 
of rotating homogeneous turbulence. The reader is referred to Cambon, Mansour & 
Squires (1994) (contained in this proceedings) for further discussion of these issues. 
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